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ABSTRACT 

We investigate the 70 /im Far-Infrared Radio Correlation (FRC) of star-forming galaxies in the 
Extended Chandra Deep Field South (ECDFS) out to z > 2. We use 70 /xm data from the Far- 
Infrared Deep Extragalactic Legacy Survey (FIDEL), which comprises the most sensitive (^0.8 mJy 
rms) and extensive far-infrared deep field observations using MIPS on the Spitzer Space Telescope, 
and 1.4 GHz radio data (~8 /iJy beam~^ rms) from the VLA. In order to quantify the evolution of 
the FRC we use both survival analysis and stacking techniques which we find give similar results. We 
also calculate the FRC using total infrared luminosity and rest-frame radio luminosity, qriR^ and find 
that qtir is constant (within 0.22) over the redshift range 0-2. We see no evidence for evolution in 
the FRC at 70 /im, which is surprising given the many factors that are expected to change this ratio 
at high redshifts. 

Subject headings: galaxies: evolution, galaxies: formation, infrared radiation, radio continuum 



1. INTRODUCTION 

The correlation between the far-infrared (FIR) and ra- 
dio emission for star-forming galaxies in the local Uni- 
verse was first observed by Ivan der KruitI (Il97ll Il97l 
and is the tightest and most universal co rrelation known 
amon g global parameters of galaxies (jHelou fc Bicavl 
119931) . The correlation is linear, spans five orders of mag- 
nitude of bolometric luminosity and has been shown to 
hold f or a wide range of Hubble types (Ide Jong et al. 

iH elou. Soifer fc R owan- RobinsonI 119851 : iCondon 
L992; Yun . Reddv fc Co ndon 20Ql|). 

The FIR-radio correlation (FRC) has been attributed 
to the presence of young, high-mass (M > 8Mo) stars. 
The FIR emission arises from the absorption by dust 
of the UV emission and subsequent reradiation of the 
energy at IR wavelengths. The radio emission is domi- 



nated by non-thermal synchrotron emission from cosmic 
ray electrons which are accelerated by supernovae shocks. 
There is also a thermal component which arises from free- 
free emission from ionized hydrogen in Hii regions, but 
this contributes only ~10% of the ra dio emission at lower 
frequencies (<5 GHz. ICondonlll992l ). and becomes more 

significant at higher frequenc ies. 

The "calorimeter theory" (|Voelkl [r989t ) suggests that 
the FRC holds because galaxies are both electron 
calorimeters and UV calorimeters so the total radio and 
IR outputs remain proportional independent of vari- 
ations within the galaxy. This theory however, re- 
quires that galaxies are optically thick to UV light from 
the young high-mass st ars, and thus does no t hold for 
optically-thin galaxies. iHelou fc Bicavl (|1993l ) proposed 
the non-calorimetric "optically thin" scenario involving a 
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correlation between disk scale height an d the escape scale 
length for cosmic ray electrons, while iBelll (|2003l ) con- 
cludes that the linearity of the FRC is a conspiracy as the 
star-formation rate is underestimated for low-luminosity 
galaxies at both radio and infrared frequencies. How- 
ever, these models typical ly leave out proton losses and 
non-s ynchrotron cooling ()Lacki. Thompson fc QuataertI 
|201C1| ). Ultimately, the physical origin of the FRC is still 
not clear. 

The far-reaching nature of the FRC has made it a 
valuable diagnostic. Some examples of its application 
include: usin g the FRC to ide ntify radio-loud AGN 
(jPonlev et al.l l2005t iNoTris et al.ll2006D : usi ng the FRC 
to define the radio luminosity/ SFR relation (|Bellll2003D : 
and, at higher redshifts, using the FRC to estimate dis- 
tances to sub millimeter galaxies without optical coun- 
terparts (e.g.. lCarilli fc Yunlll999[ ). Consequently it is of 
great importance to determine whether the FRC holds 
at high redshifts. 

The FRC may fail at high redshifts for a num- 
ber of reasons. Electrons are expected to lose energy 
by inverse Compton interactions with the cosmic mi- 
crowave background, whose energy density scales as 
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implying a lower level of radio emission at 



higher redshifts. Moreover, synchrotron emission is pro- 
portional to the magnetic field strength squared, so 
evolution of magnetic field st rength should affect the 
FRC at higher redshifts fe.g.. iMurph^ I2009D . Changes 
in the spectral energy distributions (SEDs) may also 
be expected due to evolution in dust properties and 
metallicity (e.g. . |Ainb lard et al. 2010; Hwang et al. 2010; 
ICh^ man et all l2010i) . Nonetheless, current studies 
show no firm evidence for evolution in the FRC (e.g. , 
Appleton et al.ll2004l: iSevmour et al.ll2009l: 
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Bourne et al. 2010; Ivison et al. r i2010aUbl : iSaTg^r et al.l 
IblOa b; Huynh et al. 2010). 

Herschel was launched in May 2009 and can probe 
the FIR to submillimeter regime from 55 nm to 671 /im, 
deeper than ever b e fore (iPilbratt et al.ll201(]). Mo st re- 
cently, [jirviliFin (l20i^~andl^^ used 
data from Herschel and found no evidence for evolution 
in the FRC out to z = 0.5 and z = 2 respectively. 

This paper studies the dependence of the FRC on 
redshift using deep 70 /.tm data from the Spitzer Space 
Telescope and 1.4 GHz data from the Very Large Ar- 
ray (VLA). This work differs from previous studies as 
we are using FIDEL (Far-Infrared Deep Extragalactic 
Legacy Survey) 70 fim data, which is the deepest 70 /im 
data taken to date. FIDEL reaches a point source rms 
sensitivity of 0.8 mJy at 70 /im, making it far more sen- 
sitive than pr e yious s tudies of the FIR at 70 fim such as 
I Sargent et al.l ()2010bf) whose data reached a point source 
rms sensitivity of 1.7 mJy. 

We focus on the 70 /xm data as this band probes closer 
to the dust emission peak (~100 /im) than, for example, 
24 /im. Furthermore, for 2; < 3, the 70 /xm band is not 
affected by emission from polycycli c aromatic hy d rocar - 
bons (PAHs; 7 - 12 /xm). While iBourne et all ()2010f ) 
also study the FRC using FIDEL data, their work was 
entirely based on stacking analysis. This work is the first 
to use such deep 70 /im data to study the evolution of the 
FRC based on individual sources. 

The data are described in Section 2, Section 3 describes 
the data analysis while Section 4 presents our results and 



analysis. This paper uses Hq = 71 km s ^ Mpc ^, JIm 
= 0.27 and = 0.73. 

2. DATA 
2.1. FIDEL 

FIDEL, the Far-Infrared Deep Extragalactic Legacy 
Survey, is a legacy science program (PI: Dickinson) which 
comprises the most sensitive and extensive FIR deep field 
observations using the Multiband Imaging Photometer 
for SIRTF (MIPS) on the Spitzer Space Telescope. Char- 
acteristics ^o£_the FIDE L data are described in detail 
in Magncl ireTall (|20 0^ and other papers. FIDEL ob- 
served three fields: ECDFS, ECS and GOODS-North. 
Observations were taken at 3 bands: 24 /im, 70 /im and 
160 /im, focusing specifically on the 70 /im data. This 
paper concentrates only on the 30' x 30' ECDFS field 
centered on 03 32 00, -27 48 00 (J2000). The observa- 
tions were designed to achieve roughly uniform sensitiv- 
ity at 70 /im across most of the ECDFS, although the 
FIDEL data include a somewhat deeper central region 
with data from a GO program (PI: Frayer). The mean 
70 /xm exposure time over the ECDFS is approximately 
6600s, yielding an RMS point source sensitivity of ap- 
proximately 0.8 mJy. The 24 /im exposure time varies 
considerably more over the field, from 5000 to 35000 sec- 
onds over most of the ECDFS, with an average exposure 
time of approximately 16000 seconds. The RMS point 
source sensitivity at 24 /im thus also varies, but is typ- 
ically in the range 8 to 14 /iJy over most of the field. 
FIDEL data were processed using the Mosaicking and 
Point -source Extraction (MOPEX, iMakovoz fc Marleaul 
l2005f ) package to form the mosaicked images. The final 
70 /im mosaic has a pixel scale of 4.0 "/pixel and a point 
response function (PRE) with an 18" FWHM. The final 
24 /im mosaic has a pixel scale of 1.2 "/pixel and a PRE 
with a 5.9" FWHM. 

2.2. Ancillary data 

2.2.1. Radio data 

The ECDFS has been observe d at 1.4 GHz by 
both the VLA and t he ATCA jiNorris et al.l 120061 : 
iKellermann et al.ll2008t iMiller et al.ll2008[ ). Here we use 
the lMiller et al.l Hool) radio data due to its high angular 
resolution and sensitivity over our field. The radio data 
encompass a 34' x 34' region, centered on 03 32 28.0, -27 
48 30.0 (J2000). The data have a typical rms sensitivity 
of 8 /iJy per 2.8" x 1.6" beam. The catalog contains 464 
sources above a 7 sigma cutoff. 

2.2.2. Redshift data 

We obtained both spectros copic and photom etric red- 
shift data f rom COM BO-17 (IWolf et alll2004) . MUSYC 
Gawiser et al. 200^; iCardamone et al.l I2010D GOODS 
Balestra et al.li20im and ATLAS (|Mao et aLll2009l) . 



COMBO-17 has photometric data in 17 passbands 
from 350 nm to 930 nm for 63501 objects in the ECDFS. 
The photome tric redshifts are most reliable for sources 
with R < 24 (jWolf et al.ll2T)0l . 

MUSYC (Multiwavelength Survey by Yale -Chile) has 

fotometric data in the ECDFS. Cardai none et al] 
ilO') combined photometric data from the literature 
with new deep 18-medium-band photometry into a pub- 
lic catalog of '^SOOOO galaxies in ECDFS, from which 
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they computed the photometric redshifts. The photo- 
metric redshifts_axe_jn£sti;el^ for sources with R < 
25.5. ICardamone et al.l ()2010t) also compile a spectro- 
scopic redshif t catalog of 2551 gal a xies from the the lit- 
erature (e.R., iBalestra et all 120101 : IVanzella et 111120081 : 
iLe Fevre et 311120041). 

IBalestra etal.l (|2010D observed the GOODS-South field 
(within the ECDFS) using VIMOS to obtain spectro- 
scopic redshifts. Their campaign used two different 
grisms to cover different redshift ranges and used 20 VI- 
MOS masks. They combined their resulting redshifts 
with those available in the literature to produce a cata- 
log containing 7332 spectroscopic redshifts. Quality flags 
were provided for all the redshifts and we took only those 
redshifts that had quality flags of "secure" and "likely" 
yieldin g a catalog of 5528 spec troscopic redshifts. Al- 
thoug h TCardamone et al.l ()2010D includes IBalestra et al.l 
(|2010D data in their compilation of spectroscopic red- 
s hifts, they only inc lude a subset of the data. 

iMao et al.l ([20091) are undertaking a program of red- 
shift determination and source classification of all AT- 
LAS (Austraha Telesc ope Large Area Survey) radio 
sources with AAOmega (jSharp et al.ll2006l ) on the Anglo- 
Australian Telescope (AAT). Using redshifts from both 
the literature and their own campaign they have a total of 
261 spectroscopic redshifts in ECDFS and its surround- 
ing region. 

2.2.3. X-ray data 

We use the 2 Ms Ch a ndra Deep Field South X-ray 
catalogs from iLuo et al.l ( 2008[) to identify AGN in our 
70 ^m catalog (Section l4.ip . This is one of the most 
sensitive X-ray surveys ever performed and detects 462 
sources in 436 arcmin^. While the X-ray data cover a 
smaller region than the FIDEL data, we use these data 
to eliminate some AGN from our sample (Section |43J. 

3. DATA ANALYSIS 
3.1. FIDEL 

3.1.1. 70^m 

The 70 /im catalog was produced with the Astronomi- 
cal Point Source Extraction (APEX) module within the 
MOPEX package. APEX subtracts the local background 
by calculating the median in a region, which we set to 
34 X 34 pixels, surrounding each pixel, and removing 
the 100 brightest pixels. Peak values with a S/N greater 
than three were fitted using the PRE. The 3-sigma cat- 
alog extracted using APEX contained 515 sources. 

3.1.2. g^/im 

A similar source detection and extraction process to 
the 70 /J,m data were used for the 24 /j,m data. However a 
four sigma cutoff was used to reduce the number of spu- 
rious sources. Visual inspection was required to remove 
spurious sources due to artifacts surrounding bright ob- 
jects. The final catalog of 24 /im sources in the ^30' x 
30' region of interest contained 5319 sources. 

3.1.3. Final catalog 

To test the completeness of the 70 /im catalog, we per- 
formed Monte-Carlo simulations. A simulated source 
was injected at a random location in the 70 /im im- 
age, and source extraction was performed using the same 



technique as for the production of the catalog. The in- 
put fiux density of the simulated source varied over the 
range of fluxes in the "real" image. This process was re- 
peated 10000 times and we tracked the simulated source 
recovery rate, which enables us to estimate the overall 
completeness over the entire image. Figure |T] presents 
the completeness as a function of flirx density. Our cata- 
log is almost 100% complete at 9 mJy, and 50% complete 
at ~2.5mJy. This high level of completeness at 2.5 mJy 
(^3 sigma) is probably due to the 70 /im image not being 
completely uniform. The 10 x 10 arcminute region in the 
center has an RMS point source sensitivity of ~0.6mJy, 
so 2.5 mJy is a >4 sigma limit for ■~11% of the image 
area. The completeness plot does not reach 100% due 
to the random positioning of the simulated source. If, 
by chance, the simulated source is injected upon a real 
source, the simulated source is not recovered as an indi- 
vidual source, and instead, the "real" source is recovered 
with a larger flux density. 

The 70 /im catalog of 515 sources obtained using APEX 
may contain confused sources, as well as spurious sources. 
We determined which sources required deblending by 
comparing it with the 24 /im catalog which has better 
resolution. First, we matched the 70 /im catalog to the 
24 /im catalog using a matching radius of 9" (the half- 
width at half maximum of the PRE at 70 /im). We then 
calculated the ratio of the flux density of the brightest 
24 /im source to the total fiux density of all 24 /im sources 
within the 9" matching radius (see Figure Those 
sources for which the ratio is one would not req uire de- 
blendi ng. Using arguments similar to those of Pop e et al] 
(|2006f ). we assume that if a source is bright at 24 /im, it 
will also be bright at 70 /im. Consequently, if the ra- 
tio was greater than 0.8 (that is, if the brightest 24 /im 
source contributed >80% of the total fiux density in the 
matching radius) , the 70 /im emission was determined to 
be from the brightest 24 /im source. If the ratio is less 
than 0.8, we considered the 70 /im source as a candidate 
for deblending. Using these criteria we determined that 
164 sources were deblend candidates. In principle, it is 
possible that for sources with a ratio of greater than 0.8, 
we are overestimating the 70 /im fiux density by up to 
20%. An overestimation of the 70 /im fiux density may 
result in a higher ERG, but this effect is much less than 
Agyo = 0.1. 

The sources were deblended using double, triple or 
even quadruple Gaussian fits, constraining the centers 
of the Gaussians to the 24 /im positions. We visually in- 
spected all the deblend candidates and a small number of 
these clearly did not require deblending as the secondary 
24 /im source was right on the edge of the 9 arcsecond ra- 
dius with which we computed the deblend criteria, and 
were subsequently discarded as deblend candidates. In 
total, 143 sources were deblended. resulting in a total 
catalog of 658 putative sources. In the course of the vi- 
sual inspection, we identified 41 spurious sources. 30 of 
these were sources that were faint (<4 a) in the 70 /im 
catalog, and did not have 24 /im counterparts, and 11 
were clearly part of the Airy ring of an adjacent, bright 
source. This resulted in a final catalog containing 617 
70 /im sources. 

Assuming Gaussian statistics, noise spikes are expected 
to produce approximately 34 spurious sources above our 
3-sigma cutoff in the 70 /im data, and approximately 7 
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spurious sources above our 4-signia cutoff in the 24 
data. The probability of any spurious 70 iim source lying 
within a 9 arcsec radius of any spurious 24 source is 
approximately 2%, and so we conclude that none of the 
617 sources is likely to be spurious. 

We note that after deblending the 70 /xm catalog we 
increased the number of faint sources, which changes the 
completeness levels. While we have improved our com- 
pleteness at the fainter flux levels, this change cannot 
easily be quantified. 

This catalog of 617 70 /im sources was matched to the 
24 jum catalog usin g a matching radius of 4". Following 
iHuynh et al.l ( 20081 ) , the matching radius was determined 
by plotting the number of candidate matches against po- 
sition offset (Figure [S]) . The matches from chance alone 
are determined from the source densities of the catalogs. 
Although Figure [3] suggests that a greater number of 
matches would be obtained using a matching radius of 
8" , we choose the more conservative 4" as our matching 
radius to minimize spurious matches. 

Given the relative sensitivities of the 70 fim and 24 /im 
data it is expected that most, if not all the 70 /im sources 
will have a 24 /um counterpart. A small number (10) of 
strong 70 //m sources did not have 24 /xm counterparts in 
the catalogs, while the 24 /im image showed faint detec- 
tions. We therefore performed aperture photometry on 
these faint sources and successfully extracted eight. The 
final catalog has only two 70 /im sources without 24 /im 
counterparts but neither of these have extreme S70/S24 
ratios. 

3.2. Radio, Optical and X-ray Counterparts 
3.2.1. Radio Counterparts 

The 70 nm catalog was matched to the radio catalog 
of iMiller et all (|2008l) . hereafter M08, using a matching 
radius of 5". The matching radius was determined in a 
similar manner to that described in the previous section. 
Where available, the 24 /im position was used because 
the higher resolution at 24 /im allows for better positional 
accuracy. This resulted in 171 radio sources matched to 
the 70 /im catalog. 

The M08 radio catalog has a 7 sigma cutoff. In order to 
increase the number of radio counterparts, and given that 
we have the additional information of the 70 /im sources, 
we extracted radio sources at the known 70 /im source 
positions that had radio detections greater than 3 sigma. 
This was done by performing Gaussian fits twice, first 
with a fixed size (2.8" x 1.6" - the beam of M08 data), 
and then without a fixed size. The relationship between 
peak and integrated flux densities was determined by: 



Sin 



'peak 



(1) 



where Smt is the total integrated fiux density. Speak is the 
peak flux density, Omaj is the semi-major axis of the free 
Gaussian, 9min is the semi-minor axis of the free Gaus- 
sian, bmaj is the semi-major axis of the fixed Gaussian, 
and bmin is the semi-minor axis of the fixed Gaussian. 

If Sint/ Speak was Icss than 1.2, the source was deter- 
mined to be unresolved and the peak flux density for the 
flxed Gaussian was determined to be the final flux den- 
sity. If the ratio was greater than 1.2, the source was 
determined to be extended, and the integrated flux was 



taken to be the final radio fiux density. This resulted in 
353 70 sources with radio counterparts. 

In summary, our final 70 /im catalog contains 617 
sources, 615 of which have a 24 /im counterpart. 353 of 
these have a radio detection greater than 3 sigma from 
iMiller et al.l (|2008l ). Figure H] shows the histogram of the 
70 nm flux densities. We also overplot the population 
of sources with and without radio counterparts. Sources 
with no radio counterparts have a median 70 /im flux 
density of 3.3 mJy while sources with radio counterparts 
have a median 70 /im flux density of 5.2 mJy. 

3.2.2. Redshift Data and Infrared Luminosities 

The final 70 /im catalog was matched to the various 
redshift catalogs using the matching radii shown in Ta- 
ble [H Due to the better resolution of the radio data we 
used radio positions where available, followed by 24 /im 
positions where available to match to the redshift cata- 
logs. 

Where possible we matched spectroscopic redshifts to 
the sources. Where more than one spectroscopic red- 
shift was available, we chose to prioritize ATLAS spec- 
troscopic redshifts over GOODS spectroscopic data as 
we have access to the ATLAS spectra. The difference in 
redshift, Az, was <0.0015 for the 12 sources in common. 

Where no spectroscopic data were available we used 
photometric redshifts. Where both COMBO-17 and 
MUSYC photometric redshifts were available we chose 
to prioritize MUSYC data over COMBO-17 because 
MUSYC 's photometric redshifts are derived from 32 
bands as opposed to COMBO-17, which uses 17. Fur- 
thermore, MUSYC has photometric data extending into 
the near-infrared ( JHK) , which improves the accuracy of 
the photometric redshift. 

In summary, 562 of the 617 (91%) 70 /im sources have 
redshift information, 206 (33%) of which are spectro- 
scopic (Figure O. 

The total infrared luminosity (8 - 1000 /im) of the 
sources in our sample was estimated by fitting the 
24 um and 70 /im flux d ensity to the SED templates 
of ICharv fc Elba3 (|200l . hereafter CEOl, while let- 
ting the templates scale in luminosity. The CEOl tem- 
plates show observed 24/im/70/im flux density ratios 
which are mo r e repr ese ntative of z ^ 1 galaxies than 
IDale fc Hel ou ("2002) or 'Lagach e. Dole fc Pugetl ([200l 
templates. MagncUi et al. (2009) found that 70 /im data 
seems to provide an estimate of hm that is nearly in- 
dependent of the SED library used, but this is partly 
because 70 /zm is close to the peak emission, and there- 
fore carries the largest fraction of FIR power. 

We integrate the best fit template over 8 - 1000 /im to 
derive the total IR luminosity (Lin) (Figure [6|). We find 
luminous IR galaxies (LIRGs, lO^L© < Lir < 10^'^Lq) 
are detected out to z 1.25, while ultraluminous IR 
galaxies (ULIRGs, L/fj > IO^^Lq) are detected out to z 
= 3. 

3.2.3. X-ray Counterparts 

Using the 2Ms Chandra data (|Luo et al.ll200l . we find 
55 of the 70 /im sources are within 2 arcseconds of an X- 
ray source. We use this data to calculate the hardness 
ratio so as to discriminate against AGN (Section 14. ip . 

4. RESULTS AND ANALYSIS 
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We define qm as 

g/fl = logio(-^^), (2) 

'Jradio 

where Sin is the observed infrared flux density at the 
specified IR wavelength (e.g. 70 /im), and Sradio is the 
observed flux density at 1.4 GHz. 

4.1. AGN identification 

We wish to study the FRC of predominantly star- 
forming galaxies and so we removed galaxies from our 
sample if they satisfied any of the following four criteria 
which indicate AGN. 

1. The source has a radio morphology that displays 
the classic double-lobed AGN morphlogy. Figure 
[7] shows the three sources that were identified and 
subsequently removed from our sample. 

2. The source has log(S70/S24) < 0.5, because 
AGNs are expecte d to have low S70/S24 ratios 
(jFraver et al.ll2006[ ). 22 sources were identified as 
AGN in this way. Figure [5] shows the log of the 
ratio of 70 /im fiux density over 24 /im flux density, 
plotted against redshift. 

3. The source has a ha rdness ratio > 0. 2 based on the 
2Ms Chandra data (|Luo et al.ll2008[) . where hard- 
ness ratio is defined as: 

j^j^ _ Shard — Ssoft 
Shard + Ssoft 

where HR is the hardness ratio. Shard is the flux 
density of the 2-8 keV band, Sso ft is the fiux den - 
sity of the 0.5 - 2.0 keV band (jRosati et al.ll2002[ ). 
14 sources were identified as AGN in this way. 

4. The source has a soft X-ray (0.5-2 keV) to R-band 
fiux densi t y ra tio of greater than one, as used by 
iLuo et all pOOS ) to classify AGN. 14 sources were 
identified as AGN in this way. 

Low (770 and (724 value s may also be used to discrim- 
inate against AGN (e.g. iMiddelberg et all 120081 ). How- 
ever, doing so could potentially bias our results against 
low values of q. Applying these criteria removed only a 
further five AGN from our sample and we found made a 
negligible difference to our results, so we do not use these 
criteria here. 

In total, 44 (^7%) sources were classified as AGN and 
subsequently removed from further analysis, leaving us 
with a final catalog of 573 70 fj,m sources. Most of the 
sources were identified as AGN based on only one of the 
above four diagnostics, highlighting the need for multi- 
ple AGN diagnostics. There were only six sources that 
were identified as AGN by both the X-ray diagnostics. 
41/44 of the sources identified as AGN had redshift in- 
formation. The median redshift of the AGNs is 0.969, 
slightly higher than the median redshift of the entire 
sample (0.655). The median Ljn of the AGNs is 5.10 
X 10^^ Lq, also higher than the median Ljj^ of the en- 
tire sample (2.07 x 10" Lq). The AGN tend to be at 
higher redshifts and higher Lm due to Malmquist bias. 



4.2. qjo 

To explore how the FRC, as probed by the observed 
70 /im band, changes with redshift, we plot gyo against 
redshift in Figure [SI Our data are plotted in black 
and lower limits of qjo from the radio non-detections 
are plotted as blue arrows. The lower limits are cal- 
culated by choosing the upper limit to the radio flux 
density for the radio non-detections to be 3 x rms. Data 
from the Spitzer Extr agalactic First Look Survey (xFLS, 
lAppleton et al.l |2004[ ) are plotted as red open circles. 
It is immediately evident that our data contain many 
more high redshift (z > 1) sources than the xFLS data. 
lAppleton et"all (|2004D found a qjo value of 2.15 ± 0.16 
for their sample. The data points that are at low q^Q 
values are likely to be AGN. 

We also overplot the expected ^70 ratio as a function o f 
redshift fro m IR SEP templa tes of iDale fc Heloul (|2002D 
(DH02) and |ChaQL& Elbaz! pOOll ) (CEOl). The CEOl 
SED templates are obtained from low-redshift galaxies 
with appropriate luminosities. We use four CEOl tem- 
plates which range in hjn from 10® (normal galaxy) to 
10^3 Lq (ULIRGs). The DH02 SED templates are de- 
rived from combining theoretical SEDs from dust emis- 
sion of individual regions within a galaxy. We use four 
DH02 templates with a= 1, 1.5, 2 and 2.5, where a rep- 
resents the relative contributions of the different dust 
emission SEDs. The overplotted qjo tracks are derived 
using the IR SED templates, while the radio flux is de- 
rived from the infrared flux assuming both thermal and 
non-thermal components and a constant value of q. 

In order to quantify the evolution of gyo, we binned 
all data with both 70 //m and radio detections in red- 
shift and determined the median qjo in each redshift 
bin. The results are shown in Table [2] as well as in the 
black points of Figure [TOl We calculate a median q^o of 
2.16 lb 0.03 out to z = 0.5, which is in agreement with 
lAppleton et al.l (I2004D. as well as other previous studies 
(e.g.. ISevmour et al.ll2009l ). Sources with no redshift in- 
formation have a mean gyg ratio similar to z ^ 1 sources. 

Figure Uni shows no evidence for evolution in the FIR- 
radio correlation because the data points are consistent 
with the overplotted qro tracks. At z > 1.25 we are de- 
tecting only ULIRGs (see Figure [6]), and this is reflected 
in the higher median qjo values, which arc better traced 
by the qjo track derived from the ULIRG SED template. 

4.3. Survival Analysis 

Our 70 /im selected sample has sources with both 
70 Hm and 1.4 GHz detections, as well as sources that do 
not have 1.4 GHz detections. To study the full IR sam- 
ple, we have to include radio non-detections. As we know 
the rms of the 1.4 GHz data we can estimate an upper 
limit to the radio flux density, which translates to a lower 
limit for q^Q. We choose the upper limit to the radio flux 
density for the radio non-detections to be 3 x rms. One 
method of including limits from non-detections is using 
a branch of statistics called survival analysis. Survival 
analysis is an extensive field of statistics and was first 
applied to astronomy bv lFeigelson &: NelsonI (|1985f) . We 
perform survival analysis on our data using the ASURV 
(Astronomical SURVival an alysis) package developed by 
iLavallev. Isobe fc Feigelsm] (|1992f ). 

We took the 50'^ percentile of the Kaplan-Meier esti- 
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mator to obtain a median q^o value for each redshift bin. 
The results are shown in Table [3l as well as in Figure 
[TUl The two highest redshift bins appear to have higher 
median qyo values, but this is because at z > 1.25 we are 
detecting only ULIRGs, and hence these two points lie 
along the most actively star-forming, or ULIRG track. 

Our qjQ values appear to agree, within the errors, to 
the g7o tracks derived from the empirical SED templates 
from the local Universe. This is evidence for little, if 
any evolution in gvo, which implies that galaxies at high 
redshifts (z ^ 2.5) share many of the same properties as 
galaxies in the local Universe. 

4.4. Stacking 

To extend the results we performed a stacking analysis 
of the IR sources by stacking the radio sources in bins 
of redshift. We started by stacking the M08 radio data 
at the positions of all sources not detected in the radio 
image, using a noise-weighted mean, for which the noise 
was determined from the radio image in the vicinity of 
the IR source. The stacked radio flux density was then 
obtained from a Gaussian fit to the signal at the center 
of the radio stack, leaving all Gaussian parameters un- 
constrained. The resulting stacks are shown in Figure 
[TTl and each has a significant detection. 

To derive an average qjo for all IR sources in the var- 
ious redshift bins (green points in Figure fTO|) . we calcu- 
lated a mean by adding the flux densities of the detected 
sources to the stacked flux density. Each was added with 
a weight =l/n (where n is the total number of measure- 
ments, including both stacked and detected images) and 
the stacked flux was weighted by m/n where m is the 
number of stacked images. The results are summarized 
in Table m The higher redshift bins have higher q values 
when compared to the survival analysis results because 
stacking analysis imposes lower limits on the radio flux 
densities. As before, the FIR-radio flux density ratios 
appear to follow the ratios expected for local SEDs and 
hence we find little evolution in the FIR-radio correlation 
as a function of time. 

4.5. qriR 

We computed the ratio of total infrared luminosity over 
rest frame radio luminosity to confirm our results from 
Sections 14.21 to l44l Furthermore, this allows us to make 
direct comparisons with other studies that compute the 
FRC in this way. 

We computed qriR for all 521 sources that had redshift 
information (and were not identified as AGN) using: 



QTIR = log( 



L 



IR 



3.75 X 1012W 



) - log( 



AGHz - 



WHz 



-1 



(4) 



where Ljr^ is the total infrared luminosity (Section l3.2.2|) . 
and LiAGHz is the rest-frame 1.4 GHz luminosity. In 
order to account for k-correction in the radio, we set the 
spectr al in dex, a, to 0.8^, f ollowin g the work of llbar et al.l 
((20T0I) and lSargent etall (l2010bD . 

Figure [12] shows the distribution of qriR as a func- 
tion of redshift. For sources without a radio detection 
we set an upper limit to the radio fiux of 3 x rms. 
Similar to Figure [Q] sources without a radio detection 



(ie. lower limits) tend to have higher qriR when com- 
pared to sources with a radio detection. In order to take 
into account the effect of the radio non-detections, we 
bin our data in redshift and perform survival analysis. 
Figure [13] shows the median qriR, for only the detected 
sources (black circles), and for all sources by using sur- 
vival analysis (red triangles). We have included results 
from previous studies in Figure [13] The blue starat 
z = is at qriR = 2.64, the mean qriR found bv iBelll 
((2OOI . The open upside-d own triangles are qri R. derived 
from stacking analyses by iBourne et al.l (|2010fl and the 
open squares a re the "uncorrected" qr i r, for star-forming 
galaxies, from ISargent et al.l (|2010bl ) (See Section 14.5.11 
for a discussion on the Kellermann Correction). 

Our qriR values are all within '^0.22 of each other. 
Our values agree, wi thin the errors, with the work of 
Sargent et al.l (l2010bD . and also agree with the work of 



Bourne et al.l (j201Cl[ ). with the exception of the redshift 



1 S. 



where is the flux density at frequency u. 



bins 0.75 < z <1 and z > 1.5, where ou r values differ 
by <2 cr compared to the values gi ven bv [Bou rne et al.l 
(|2010f ) for similar redshift ranges. iSargent et al.. (2010bf) 
use 3 X rms to determine the upper limit of the ra- 
dio flux dens i ty for the ir dataset. Unlike th e work by 
IBourne et all (|2010[ ) and lSargent et al.l (|2010bD , our qriR 
values do not consistently increase or decrease over the 
redshift range studied, although we do see marginally 
significant evidence for a decreasing qriR. over the red- 
shift range to 1.5, followed by an increase from 1.5 
to 3. However, Fig 14 shows that the data points that 
use survival analysis (which are more reliable than those 
that do not) are statistically consistent with a horizon- 
tal line at about qriR = 2.62, and so we conclude that 
these apparent variations of qTiR with redshift may not 
b e statistically significant. 

'Bourne e t al.l (j2010) find their qriR values are sys- 
tematically higher than the median qriR, found bvlBelil 
(|2003( 1 of 2.64 ± 0.02, except at z > 1 where it appears 
to decline. They suggest the apparent decline may be 
attribut ed to the assumptions they make about spectral 
indices. ISargent et al.l ((2010b') find that their qxiR values 
are constant with redshift when they apply a correction 
(see Section 14.5.111 , otherw i se ~0 .3 dex of evolution is 
found. While B ourne et al.l (|2010( ) also uses FIDEL data 
in the ECDFS, their sample is IRAQ selec t ed whe reas 
our sample is 70 /xm selected. ISargent et al.l (|2010b[ ) use 
a 24 /im selected sample to study the qriR in the COS- 
MOS field. 

There are a number of factors that may contribute to 
uncertainties in qTiR- The predominant source of uncer- 
tainty is in the estimation of Lir as we are fitting the 
24 /im and 70 /zm fluxes to SED templates. The uncer- 
taint y from the choice of model SEDs can add .2 to 0.3 
dex iMagnein et al.ll2009l: [LeFloc'h et al.ll2005D . Uncer- 
tainties also arise from the rest-frame radio luminosity 
w here we are assumi ng a = 0.8. 

iJarvis et al.l (|2010( ) suggest that the slight upturn in 
qiR seen at low r edshifts (z < . 5) i n high redshift studies 
of the FRC (e.g. IBourne et al.l[201Q) may be attributed 
to the resolving of extended structure, which means ra- 
dio flux would be missed and hence increase the FRC. 
Our radio data has a 2.8" x 1.6" beam, so if the Jarvis 
hypothesis is correct, we would expect to see this effect 
at low redshifts. 

Although we see a slight upturn in qxiR from z = 1 
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to z = 0, this increase has a roughly constant gradient, 
implying that resolution effects are important as high as 
z = 1. However, at z = 1 a typical star- forming galaxy 
of diameter 10 kpc is unresolved by our beam. Thus, the 
increasing qtir towards low redshifts is probably not due 
to this effect. 

4.5.1. The Kellermann Correction 

iKellermannI ()1964l) showed that flux ratios (such as q) 
or spectral indices of a flux- limited sample are biased by a 
factor which depends sensitively on the scatter in the flux 
ratio and on the source intensity distribution, and in par- 
ticular on the power law index (3 of the differential source 
counts (i.e., dN/dS oc ). Other formulations of the 
same effect are given bv lCondonI (flQSl . IfVancisl ([l99l 
and iLauer et al.l ()2007t ). and w ere first noted as bein g 
r elevant to the evolution of q bv lSargent et al.l ()2010aD . 

ISargent et all (12010b') assume a Euclidean (/3=2.5) 
source intensity distribution at z < 1.4 and a sub- 
Euclidean (/3 ^1.5) at 1.4 < z < 2 to derive a correc- 
tion of 0.22 in their value of q at z > 1.4. However we 
note that (a) the expression assumes a common value 
of P for both the radio and infrared source counts, and 
(b) the value of /3 is very uncertain at low flux densities. 
For example, the median radio flux density for our de- 
tected sources at z < 1.4 is 72 /zJy, at which the radio 
sourc e count power law i ndex lies in the range ^1.5-2.5 
(e.g., iHuvnh et al.ll2005l ). leading to a correction factor 
between 0.22 (the adopted value of Sargent, which gives 
no evolution of q with z) and 0.51 (which implies that q 
decreases with z with marginal significance). 

Because of this uncertainty in the value of the Keller- 
mann correction for faint sources such as those stud- 
ied here and for consistency with other authors (e.g., 
Appleton et al.ll2004t ISevmour et al"1l2009t iBourne et alJ 
2010( ). we do not apply it to our data but note that it 
is responsible for a further uncertainty in the slope of q 
as a function of redshift. The range of possible Keller- 
mann corrections (which depend strongly on the assumed 
value of /3) to the variation of q over our redshift range 
is roughly centered on zero, with a total spread of about 
± 0.3. 

5. SUMMARY AND CONCLUSIONS 

We have studied the FRC out to z > 2 of ULIRGs in 
ECDFS. Our results for q-jo showed that they could be 
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broadly desc ribed by q-^q tra c ks de rived from the SEP 
templ ates of ICharv fc Elba3 (|2001l ) and iDale fc Heloiil 
(|2002| ). To quantify the evolution of q^o we binned our 
data in redshift and determined the median qyo using 
both survival analysis and a stacking analysis. Both sur- 
vival analysis and the stacking analysis gave similar re- 
sults. We see no clear evidence for evolution in the FRC 
at 70 /im. 

We also calculate the FRC using Ljji and Li ^oHz and 
find that evolution in qTiR is constrained within 0.22. 
Our calculated qriR appears slightly lower than previ- 
ous studies but this may merely be due to uncertainties 
involved in calculating qriR- We also acknowledge the 
importance of the Kellermann correction but due to un- 
certainties in the value we do not apply it to our data. 

A lack of evolution in the FRC is surprising because it 
implies that the myriad of effects on which the FRC rely 
must either also not evolve, or must evolve in such a way 
so as to preserve the FRC. E ffects such as i nverse Comp- 
ton cooling of the electrons (lMurphvll2009f ). evolution of 
the magnetic field strength, and evolution in the SEDs 
due to dust and metallicity changes should all affect the 
FRC. 

We therefore conclude that either all these factors are 
insignificant at z ~ 2 or there is a complex interplay 
between these factors conspiring in the preservation of 
the FRC at high redshifts. 

Early science results from Herschel have already hinted 
that the FRC shows no evidence for e volution to z = 2 
(jJarvis et al.ll2010H!vison et al.ll2010bl ). More results on 
the FRC from Herschel are expected over the next few 
years as observations are completed and the data ana- 
lyzed. Herschel will measure the far-infrared properties 
of normal galaxies to z 1 and ULIRGs out to z ~ 4. 
Consequently, we will be able to study the FRC out to 
higher redshifts and hence gain a better understanding 
of the evolution of star- forming galaxies, especially in the 
high redshift Universe. 
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TABLE 1 

Table of redshift catalogs. Column 1 names the catalog while Column 2 gives the matching radius which was used to 
match the redshift data to the 70 /.tm data. columns 3 and 4 give the catalog size and whether the redshift data is 
photometric or spectroscopic and columns 5 and 6 give the total number of sources that are matched to the 70 /.tm data 
as well as the number that are in the final catalog. column 7 gives the reference for the catalog. there are 562 

sources with redshifts, of which 206 are spectroscopic. 



Catalog Radius (arcsec) Spec/Phot Size ^matched ^ final Reference 



COMBO-17 


2 


phot 


62337 


555 


122 


Wolf et al. (2004) 


MUSYC 


1.5 


phot 


59693 


410 


234 


Cardamone et al. (2010) 


MUSYC 


1.5 


spec 


2551 


114 


26 


Cardamone et al. (2010") 


GOODS 


2 


spec 


7332 (5528) 


160 


148 


Balcstra ct al. (2010) 


ATLAS 


5 


spec 


254 


32 


32 


_Mao ct al. (2009) 



TABLE 2 

Median qjQ for different redshift bins for only sources with both 70um and radio detections, and are not classed as AGN. 

There are a total of 323 sources. 



z 


median z 


sources 


median qyo 


< 2 < 0.25 


0.148 


53 


2.20 ± 0.05 


0.25 < 2 < 0.50 


0.370 


46 


2.07 ± 0.04 


0.50 < 2 < 0.75 


0.626 


93 


1.87 ± 0.04 


0.75 < 2 < 1.00 


0.890 


43 


1.81 ± 0.03 


1.00 < 2 < 1.50 


1.13 


52 


1.80 ± 0.04 


2 > 1.50 


1.91 


10 


1.65 ± 0.11 


no redshift 




26 


1.79 ± 0.06 



TABLE 3 

Median gro for different redshift bins, calculated using survival analysis, for all 70um sources that are not classed as 
AGN. Survival analysis was performed for radio flux density limits of 3 x rms. The range of values given for the gyo 
values derived from model bed templates of DH02 ARE for a = 2.5 (first value) and a = 1. The range of values given for 

the model (770 VALUES DERIVED FROM SED TEMPLATES OF CEOl ARE FOR NORMAL GALAXIES AND ULIRGS. THERE ARE A TOTAL OF 573 

SOURCES. 



2 median z Nsources median 1770 (3(7) qro (DH02) qro (CEOl) 



< 2 < 0.25 


0.141 


107 


2.39 ± 0.06 


2.11 - 


2.32 


1.87 - 


2.22 


0.25 < 2 < 0.50 


0.370 


80 


2.20 ± 0.06 


2.00 - 


2.35 


1.90 - 


2.13 


0.50 < 2 < 0.75 


0.646 


154 


2.06 ± 0.04 


1.88 - 


2.33 


1.88 - 


2.01 


0.75 < 2 < 1.00 


0.862 


75 


2.01 ± 0.06 


1.77 - 


2.28 


1.83 - 


1.92 


1.00 < 2 < 1.50 


1.13 


85 


1.97 ± 0.06 


1.66 - 


2.20 


1.75 - 


1.81 


2 > 1.50 


1.88 


20 


1.90 ± 0.12 


1.41 - 


1.92 


1.48 - 


1.55 


no redshift 




52 


2.18 ± 0.08 











TABLE 4 

Summary of radio stacking results of IR sources. N^et is the number of IR sources in each bin which are detected in the 

RADIO IMAGE. Ngtack IS THE NUMBER OF IR SOURCES NOT DETECTED IN THE RADIO IMAGE AND THEREFORE THE NUMBER STACKED IN THIS 
analysis, RESULTING IN THE STACKED SuGHz LISTED IN COLUMN 5. AVERAGE (J70 IS THE FIR-RADIO FLUX DENSITY RATIO FOR ALL IR 
SOURCES IN THE REDSHIFT BIN, OBTAINED BY COMBINING THE DETECTIONS WITH THE STACKED RESULTS. 



2 median 2 Ng^t N'stack stacked Sia GHz (fJ-^j) average (J7ci 



< 2 < 0.25 


0.141 


53 


54 


27.7 ± 2.8 


9 1 1 +u.(K 
^■^^-0.03 


0.25 < 2 < 0.50 


0.370 


46 


34 


25.4 ± 1.9 


9 19+0.03 
''•^''-0.04 


0.50 < 2 < 0.75 


0.646 


93 


61 


19.7 ± 1.5 


1 at:+0.04 
^■''■:'-0.04 


0.75 < 2 < 1.00 


0.862 


43 


32 


18.9 ± 1.6 


^■»°-0.06 


1.00 < 2 < 1.50 


1.13 


52 


33 


13.6 ± 0.9 


^■**y-0.05 


2 > 1.50 


1.88 


10 


10 


19.5 ± 2.5 


1 71+0.04 
' J^-0.04 


no redshift 




26 


26 


46.8 ± 9.4 


1 Qo+0.06 
^■''■^'-O.OT 
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TABLE 5 

Median qtir for different redshift bins for only sources with both 70 /im and radio detections, and are not classed as 

AGN. There are a total of 297 sources. 



2 median z Ngaurces median Qtir 



<z < 0.25 


0.148 


53 


2.54 ± 0.05 


0.25 < 2 < 0.50 


0.370 


46 


2.52 ± 0.04 


0.50 <z< 0.75 


0.626 


93 


2.41 ± 0.03 


0.75 <z< 1.00 


0.890 


43 


2.35 ± 0.04 


1.00 < 2 < 1.50 


l.l.'^ 


52 


2.33 ± 0.03 


z > 1.50 


1.91 


10 


2.52 ± 0.11 



TABLE 6 

Median qtir for different redshift bins, calculated usinc; survival analysis, for all 70 /im sources that are not classed as 
AGN. Survival analysis was performed using Li,4ghz calculated from radio flux density limits of 3 x rms. There are a 



total of 521 sources. 


z 


median z 


^ sources 


median qtir (3(t) 


0<z< 0.25 


0.141 


107 


2.74 ± 0.06 


0.25 <z< 0.50 


0.370 


80 


2.70 ± 0.06 


0.50 <z< 0.75 


0.646 


154 


2.60 ± 0.04 


0.75 < 2 < 1.00 


0.862 


75 


2.53 ± 0.06 


1.00 < 2 < 1.50 


1.13 


85 


2.55 ± 0.06 


2 > 1.50 


1.88 


20 


2.75 ± 0.12 
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Fig. 2. — Histogram of the ratio of the maximum 24 /im flux density over the total 24 /xm flux density within a 9" radius. Sources with 
a ratio less than 0.8 were deemed deblend candidates. 
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Fig. 4. — Histogram of 70 fim fluxes. The black line is the histogram for all 617 sources, while the red diagonally shaded histogram shows 
the 353 sources that have a radio detection and the blue horizontally shaded histogram shows the 264 sources that do not have a radio 
detection. 
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Fig. 5. — Histogram of all 562 sources that have redshift information. The black line is the histogram for all sources that have redshift 
information while the red diagonally shaded histogram shows sources with spectroscopic redshifts and the blue horizontally shaded histogram 
shows sources with photometric redshifts. 
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Fig. 6. — The total IR luminosity (8 - 1000 fim) as a function of redshift for all IR sources with redshift information (562/617). Sources 
with spectroscopic redshifts are shown with filled circles while sources with photometric redshifts are shown with open circles. 
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Fig. 7. — Sources which wo classified as AGN based on radio morphology. The greyscale image is the 70 ^m image and contours are the 
1.4 GHz image starting at 24 /ijy (~3 times the rms) and increasing by factors of 2. 




Fig. 8. — The ratio of 70 fim flux density over 24 fim fiux density plotted against redshift. Sources with no redshift information are shown 
at an artificial redshift of -0.05. The blue line represents a simple modified blackbody BED model for starbursts with a dust temperature 
of 30K with a mid-infra red slope of a = 2 .4, while the red line represents a model for AGNs with a dust temperature of 90K with a 
mid-infrared slope of 1.1 lIFraver et al.ll200^ . Error bars are derived by combining in quadrature the standard errors in radio and IR fluxes. 




Fig. 9. — The 70 fira FIR-radio correlation plotted against redshift. The black points are sources that have both 70 fim and 1.4 GHz 
detections while the b lue arrows show the lower limit of qro for sources with no radio detection. The red open circles are data from xFLS 
l|ApDleton et al.ll200 4). The black lines are expected qro tracks derived from SED templates for galaxies with total infrared luminosities of 
10^ (normal galaxies), lO^^, lO^^ and lO^^Lp, (ULIRGs) going from bottom to top, from lCharv fc Elbaj 11200 II ) and the blue lines are qyo 
tracks derived from SED templates from Dale & Hclou (2002). Sources with no redshift information are shown at an artificial redshift of 
-0.05. Error bars are derived by combining in quadrature the standard errors in radio and IR fluxes. 
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Fig. 10. — Median q values for different redshift bins. The black circles show the median qro values for sources with both 70 nm and 
1.4 GHz detections, the red triangles show the median qro value for all sources taking into account the lower limits using survival analysis 
and the green squares show the median qjo value after stacking. The tracks derived from the SED templates are the same as in Figure [9] 
Vertical error bars denote standard errors and horizontal error bars indicate the range of the redshift bin. 
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Fig. 11. — Postage stamps of the radio stacks of the IR sources, for the redshift bins as shown. The stacks are 16 X 16 arcsec in size. 
There is significant radio flux at the center of the stack in all of the redshift bins. The measured flux densities arc reported in Table (4] 
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Fig. 13. — Median qtir values for different redshift bins. The black filled circles show the median qxiR values for sources with both 
70 fim and 1.4 GHz detections and the red triangles show the median qtir value for all sources t aking into account th e lower limits using 
survival analysis. The grey upside-down open tria ngles show the median qtir values derived by Bour ne et al.l l|201CII ') and the grey open 
squares show the median Qti r valu es derived by ISargent et al.l l l201 0b) for star-forming galaxies. None of the data shown include the 
Kellermann correction ('Section l4.5.H l. Vertical error bars are standard errors for both our dataset and Bourne's data, but the vertical error 
bars for Sargent 's data are upper and lower 95% confidence levels. The blue star at (l+z) = 1 represents the median qTiR= 2.64 ± 0.02 
from lBelil | |200^ . 



